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Abstract We consider an open one dimensional lattice gas on sites % = 
1, . . . ,N, with particles jumping independently with rate 1 to neighboring 
interior empty sites, the simple symmetric exclusion process. The particle 
fluxes at the left and right boundaries, corresponding to exchanges with 
reservoirs at different chemical potentials, create a stationary nonequilibrium 
state (SNS) with a steady flux of particles through the system. The mean 
density profile in this state, which is linear, describes the typical behavior 
of a macroscopic system, i.e., this profile occurs with probability 1 when 
N — > oo. The probability of microscopic configurations corresponding to 
some other profile p(x), x = i/N, has the asymptotic form exp[— NF({p})]] 
T is the large deviation functional. In contrast to equilibrium systems, for 
which J~eq({p}) is just the integral of the appropriately normalized local 
free energy density, the J- we find here for the nonequilibrium system is a 
nonlocal function of p. This gives rise to the long range correlations in the 
SNS predicted by fluctuating hydrodynamics and suggests similar non-local 
behavior of T in general SNS, where the long range correlations have been 
observed experimentally. 
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1 Introduction 



The extension of the central object of equilibrium statistical mechanics, en- 
tropy or free energy, to stationary nonequilibrium systems in which there is 
a transport of matter or energy has proved difficult. One knows from vari- 
ous approximate theories like fluctuating hydrodynamics that such systems 
exhibit long range correlations very different from those of equilibrium sys- 
tems jy, g, g]. These correlations extend over macroscopic distances, as has 
been established rigorously in some models, and reflect the intrinsic non- 
additivity of such systems. They have been measured experimentally in a 
fluid with a steady heat current j|, Q . Their derivation from a well defined 
macroscopic functional valid beyond local equilibrium (where there are no 
such correlations) is clearly desirable. We report here what we believe is the 
first exact derivation of such a functional for a nonequilibrium model which 
is relatively simple but exhibits the realistic feature of macroscopically long 
range correlations. 

Before describing our model and results, we review briefly the corre- 
sponding results for equilibrium systems || ^, 0, §. Let us ask for the 
probability of finding an isolated macroscopic equilibrium system, having a 
given number of particles and energy and contained in a given volume V, in 
a specified macro state M, that is, of finding the microscopic configuration 
X of the system in a certain region Tj^ of its phase space. According to the 
basic tenet of equilibrium statistical mechanics, embedded in the Boltzmann- 
Gibbs-Einstein formalism, this probability is proportional to exp[5'(A / i)/fcs], 
where the (Boltzmann) entropy S(A4) of the macro state M is defined by 
S(M) = ks log |r^vf | , with \Tm \ the phase space volume of T_m || ^, [7|, ||, 0]. 

When the system is not isolated but is part of a much larger system — 
a situation idealized by considering the system to be in contact with an 
infinite thermal reservoir at temperature T and chemical potential v — then 
the entropy in the formula above is replaced by an appropriate free energy. 
Consider in particular a lattice gas in a unit cube containing L d sites with 
spacing 1/L (similar formulas will hold for continuum systems), and suppose 
that the macro state of interest is specified by a density profile prescribing 
the density p[x) at each macroscopic position x in the cube. Then the 
probability of finding the system in this macro state is given for large L by 



P({p{x)}) ~ exp [ - L d F eq ({p})} 



(1.1) 



with 




(1.2) 
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The integration in ( |1.2| ) is over the unit cube, f v (f) = o( r ) — vr, where a(r) 
is the usual Helmholtz free energy density for a uniform equilibrium system 
at density r, and the equilibrium density p = p(y) corresponding to the 
chemical potential v is obtained by minimizing f u : 

da(r) 



dr 



(1.3) 

--p(u) 



Note that —f v {p) is just the pressure in the grand canonical ensemble. We 
have suppressed the dependence of / and JF eq on the constant temperature 
T, and assume for simplicity that neither p nor p(x), for any x, lies in a 
phase transition region at this temperature. 

The challenge is to extend these results to nonequilibrium systems, in 
particular, to systems which are maintained in a stationary nonequilibrium 
state (SNS) with a steady flux of particles by contact with two boundary 
reservoirs at different chemical potentials. We would like to generalize the 
formula (1.2) for JF oq ({p}), obtaining a large deviation functional J-({p}) 
such that the probability of observing a density profile p{x) is given by a 
formula analogous to (|P|). The typical profile would then correspond to 
the p which minimizes T . A similar analysis and an appropriate J- would 
certainly be useful for the study of pattern formation in more general SNS. 
An interesting example is the Benard system, in which the particle flux is 
replaced by a heat flux maintained by reservoirs at different temperatures, 
with the hotter reservoir below the system, and the typical patterns change 
abruptly from uniform to rolls to hexagonal cells as the fluid is driven away 
from equilibrium [jl(|. Open systems of this type have been discussed ex- 
tensively in the literature from both a microscopic and macroscopic point of 
view; see [|ll], 12] and references therein. 



Given our current limited understanding of such SNS, however, it is 
necessary to start with the simplest systems; our results here are for the 
one dimensional symmetric simple exclusion process (SSEP) [^, [14|, p"5| ], 
driven by boundary reservoirs at distinct chemical potentials vq and v\. We 
consider a lattice of N sites, in which each site i is either empty (t{ = 0) 
or occupied by a single particle (tj = 1), so that each of the 2^ possible 
configurations of the system is characterized by N binary variables ri, rjy- 
Each particle independently attempts to jump to its right neighboring site, 
and to its left neighboring site, in each case at rate 1 (so that there is 
no preferred direction). It succeeds if the target site is empty; otherwise 
nothing happens. At the boundary sites, 1 and N, particles are added or 
removed: a particle is added to site 1, when the site is empty, at rate a, 



3 



and removed, when the site is occupied, at rate 7; similarly particles are 
added to site ./V at rate 5 and removed at rate (3. This corresponds to 
the system being in contact with infinite left and right reservoirs having 
fugacities zq = expz^o = a/7, and z\ = expi^i = 5/(3, see [[[(], [H]]. We 
therefore define 

z a zi 5 

1 + zq 7 + a I + 21 (3 + 5 

and think of these as the densities of the reservoirs. They will in fact be the 
stationary densities at the left and right ends of the system when N — > 00 
sec O below). 

We refer to zq and z\ as fugacities because if we were to place our system 
in contact with only the left (right) reservoir by limiting particle input and 
output to just the left (right) side, i.e., by setting 5 = (3 = (a = 7 = 0), 
then its stationary state would be one of equilibrium with fugacity zq (zi), 
with no net flux of particles. Of course if zq = Z\ = z then the system would 
be in an equilibrium state whether in contact with one or both reservoirs, 
i.e., in a product measure with uniform density p = z/(l + z)\ this value 
follows from (|0]), since for this system 

a(r) = r log r + (1 — r) log(l — r). (1-5) 

For the model considered here the typical profile p(x,t), on the macro- 
scopic spatial-temporal scale with variables (x,t) defined by i — > xN and 



tN 2 , is for N —¥ 00 governed by the diffusion equation Jl6|, |17 
dp(x,t) d 2 p(x,t) 



dt dx 2 



p(0,t)=p o , p(l,t) = Pu (1.6) 



which gives for the stationary state p(x) = po + (pi — po)x. Spohn || (see 
also [|l!|) computed explicitly, for a + 7 = + 6 = 1, both the expected 
density profile and the pair correlation in the stationary state, and showed 
that the results agree, in the above scaling limit, with those obtained from 
fluctuating hydrodynamics. 

While higher order correlations can also be obtained in principle, the 
difficulty of the computation increases rapidly with the desired order. The 
complete measure on the microscopic configurations in the steady state may, 
however, be computed through the so-called matrix method [^0|, [2l], |22|, [23| , 
From this measure we would like to determine the probability of seeing 
an arbitrary macroscopic density profile p(x),x G [0,1]; by definition, this 
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is the sum of the probabilities of all microscopic configurations which are 
consistent with p(x), e.g., all configurations r such that for any xq,xi with 
< xq < xi < 1, jf J2iix Q N T i ~ Ixo P( x ) < for some appropriate 
choice of 5n with 5n — > as N — > oo, e.g., <5jv = A r_K , < re < 1. 

The deviations of macroscopic systems from typical behavior, e.g., from 
the profile p(x, t) governed by the diffusion equation, is a central issue of 
statistical mechanics || ^, 0, |8[ |ll], |l2| . The time dependent problem was 
first studied by Onsager and Machlup p5| , who considered space-time fluc- 
tuations for Hamiltonian systems in equilibrium (for which p is constant). 
The problem of observing a profile p(x,t) not necessarily close to p(x,t), 
the so-called problem of large deviations, was later studied rigorously for 
the SSEP on a torus by Kipnis, Olla, and Varadhan |2q| , who found an 
exact expression for the probability of observing such a profile (see also 
jTBL [l8|]). An important ingredient in their analysis was the fact that the 
dynamics satisfied detailed balance with respect to any product measure 
with constant density, and that such a product measure is in fact the appro- 
priate microscopic measure at each macroscopic position. This corresponds 
for stochastic systems to reversibility of the microscopic dynamics, which 
also plays an important role in the Onsager-Machlup theory. There have 
been various attempts to extend these results to open systems in which the 
dynamics no longer satisfy detailed balance with respect to the stationary 
measure [jTJ, 27, 28, |2{|. Bertini et al. pl| have indeed succeeded in doing 



this for the so-called zero-range process with open boundaries, for which the 
microscopic stationary state is also a product measure. The problem is more 
difficult for our system, where it is known |J that the stationary microscopic 
state contains long range correlations — correlations expected to be generic 
for SNS @, H, gl], H and even measured experimentally in some cases 



We shall see that this is reflected in a nonlocal structure for ^({p}). A pre- 
liminary version of our results was given in [33|, and more recently Bertini et 
al. [29] succeeded in rederiving this result, and obtaining a large deviation 
functional for space-time profiles, using a semi-macroscopic method (which 
could be further checked, at least for small fluctuations, by comparing its 
predictions with the direct calculation of the time dependent correlations 
given in [34]). 
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2 Summary of results 



In the present paper we give the exact asymptotic formula for the proba- 
bility Pn({p(x}) of seeing a density profile p(x), < p(x) < 1, in the one 
dimensional open system with SSEP internal dynamics. (The exact sense 
in which we establish this formula is sketched later in this section.) Let po 
and p\ be defined as in (|1.4| ); we will assume for definiteness that po > pi 
(results for the case po < Pi then follow from the right-left symmetry or 
the particle-hole symmetry, and for the case po = p\ by taking the limit 

Po -» Pi)- 

Our main result is then: 

hm l0gP \ ( ^ (x}) = -T({P}) , (2-1) 



where 



Hip)) = f Q dx{p(x)\og(^) + (i- P (x))\or (l p{x 



1 - F(x) 



+ (2.2) 



Pi - Po 



The auxiliary function F(x) in (^2) is given as a function of the density 



profile p(x) by the monotone solution of the nonlinear differential equation 

with the boundary conditions 

F(0) = po , F(l) = Pl . (2.4) 

We will show in Section |5| that such a solution exists and is (at least when 
po < 1 and pi > 0) unique. 



We now summarize some consequences of (pO -2.4). 
(a) Let us denote the right hand side of (|2.2| ), considered as a functional 
of two independent functions p(x) and F(x), by Q: 

p( x ) \ ( 1 — p( x 



G({P},{F})) = jf 1 cfa{p(x)log(|^) + (l-p(a J ))log 



F(x) 



log ( ) l . ( 2 .5) 



Pi ~ Po 
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If one looks for a monotone function F, satisfying the constraint (2.4), for 



which G{{p}, {F}) is an extremum, one obtains fl2.3|) as the corresponding 
Euler-Lagrange equation: 

m{p},{F}) _ 

oF(x) 

We will show in Section || that the unique monotone solution of ( |2,3| ) (or 
equivalently of ( |2,6D ) is in fact a maximizer of (|2.5| ): 



H{p})=^pG({p},{F}) . (2.7) 

We will also show that T is a convex function of {p}. (More precisely, 
our derivations will be for the case 1 > po > p\ > 0, but we expect the 
conclusions to hold in general, that is, for 1 > po > Pi > and by symmetry 
for 1 > pi > p > 0.) 

(b) From (|2.6D one sees at once that for F the solution of ( |2.3D , (|2~1 



5T({ P }) _ 5g({ P },{F}) _ bg 

<5p(x) ^p( x ) 



p(x) 1 - F(x) 
F(x) ' 1 - p(x) 



(2. 



and this together with ( |2.3| ) implies that the minimum of ^({p}) occurs for 
p(x) = F(x) = p(x), where 

p[x) = p (l - x) + pxx. (2.9) 



Moreover, from fl2.2| ) and (2.3) one has ^"({p}) = 0, confirming that the most 
likely profile p(x) is obtained with probability one in the limit N — > oo. Any 
other profile will have ^({p}) > and thus, for large N, exponentially small 
probability. The profile may be discontinuous, or may fail to satisfy the 
boundary conditions p(0) = po or p(l) = p%, and still satisfy T{{p\) < oo. 
When po = or pi = 1 there are some profiles for which T = +oo; their 
probability is super-exponentially small in N. For examples, see (d) below 
and Section |8[ 

(c) It is natural to contrast the SNS under consideration here with a 
local equilibrium Gibbs measure for the same system — a lattice gas with 
only hard core exclusion — with no reservoirs at the boundaries but with 
a spatially varying chemical potential u(x) |2^, ^] which is adjusted to 
maintain the same optimal profile p(x). For this system the large deviation 
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functional (free energy) is obtained directly from Ql.2| ), with f u {r) = r logr+ 
(1 — r) log(l — r) — vr: 



p( x ) bg ?W + [1 - p {x)\ log (1 \ ,h: 



(2.10) 



p(x) (l-p(x)) 
In general the two expressions (^^) and ( |2.10| ) are different, and from ( |2.7| ), 
^({p}) = sup£({p},{F}) > g({p},{p})=^({p}) • (2-11) 

F({p}) and ^eq({p}) agree only for p(x) = p(x) or in the limiting case 
Po = Pit i n which the system is in equilibrium with p{x) = p$. (Equations 
(|2,2|-|2.4|) have a well-defined limit for p\ f p 0; with F(x) = Po + (p± — 
po)x + 0((pi — po) 2 )-) Otherwise F({p}) lies above ^" e q({p}) and thus gives 
reduced probability for fluctuations away from the typical profile. 

Note that the integrand in ( |2.10| ) (or ( |1.2|) ) is local: changing p(x) in 
some interval [a, b] only changes the value of this integrand inside that in- 
terval. This is not true for the integrand of ( |2.2[) , because F is determined 
by the differential equation (2.3) and so F(x) will generally depend on the 
value of p(y) everywhere in [0, 1]. 

(d) For a constant profile p(x) = r, the solution F of (p.3|) and (|2.4 
satisfies F' = AF T {1 - F) l ~ r , where A is fixed by (|J), and 



H{P}) = log 



in 



I'o 



1 



1-r 



dz 



Pi ~ Po 



(2.12) 



We see that T({p\) = oo if r = and po = 1, or r 
contrast ^-" e q({p}) as given in (|2TTo|) would be 



1 and pi = 0. By 



^eq(M) 



Pi 



log 



/'(i 



1 



dz 



Pi ~ Po 



(2.13) 



which is finite except in the degenerate cases po = p\ = 1 and po = pi = 0, 
in which case the measure is concentrated on a single configuration and 
^"eq = 00 unless p = po- 

(e) If we minimize T subject to the constraint of a fixed mean density 
Jq 1 p(x)dx, the right hand side of ( |2.8D becomes an arbitrary constant, and 
together with (|2.3[) one obtains that the most likely profile is exponential: 
p(x) = A x exp(0x) + A 2 (with F(x) = 1 - A 2 + A 2 (l - A 2 ) exp{-6x)/A 1 ), 
the constants being determined by the value of the mean density and the 
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boundary conditions fl2.4)- (This exponential form, which is the stationary 
solution of a diffusion equation with drift, was first suggested to us by Errico 
Presutti). 

Similarly, if we impose a fixed mean density in k nonoverlapping inter- 
vals, p{x)dx = Cj for i = 1, . . . , k, with no other constraints, then one 
can show using (^J^) and ( |2.8D that the optimal profile has an exponential 
form inside these intervals and is linear outside; it will in general not be 
continuous at the end points of the intervals. 

(f) Using the fact that the exponential is the optimal profile for the 
case of a fixed mean density in the entire interval, we may compute the 
distribution of M, the total number of particles in the system, in the steady 
state for large N. We find that the fluctuations of M predicted b y (|2,2| ) are 
reduced in comparison to those in a system in local equilibrium ( 2.10| ) with 
the same p: 

Jim A^KM^sNS-iM} 2 ] 

= lim A" _1 [(M 2 ) eq - (M) 2 ] - iP±^J*t. (2.14) 

N—>oo LZ 



(Since p is given by (]279|) for both systems, (M)sns = (M) C(l = (M) = 
(po + Pi)/2-) This reduction of fluctuations was already visible in fl2.11p . 
We may also obtain (|2.14| ) by expanding p(x) about p(x) in (|2.2|) . The 
result agrees with that obtained in Q directly from the microscopic model 
and from fluctuating hydrodynamics 0. 

The structure of the rest of the paper is as follows: In Sections || and 
U and Appendices [A] and ^ we derive our main result fl2.2| - |2.4| ) from the 
knowledge of the weights of the microscopic configurations as given by their 
matrix product expressions. Our approach consists in considering the sys- 
tem of N sites as decomposed into n boxes o£ Ni, N2, ■ ■ ■ , N n sites. We first 
calculate the generating function of the probability PNi,...,N n (Mi, • • • iM n ), 
that Mi particles are located on the first Aq sites, M% particles on the next 
N2 sites of the lattice, etc. We then use this generating function to obtain 
PNi,...,N n (Mi, . . . , M n ). The expression we obtain appears in a paramet- 
ric form because we extract this probability from the generating function 
through a Legendre transformation. Then we take the limit of an infi- 
nite system which corresponds implicitly to the usual hydrodynamic scaling 
limit, as explained in great detail in flq, by first letting — > 00, keeping 
Ni/N = yi fixed, and then letting j/j — > to obtain (^2-2.4). 



In Section |5| we prove that for any profile p{x) there exists a unique 
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monotonic function F(x) which satisfies ( |2.3| ) and Q2.4| ). We also establish 
there, for 1 > po > p\ > 0, that T{{p}) = sup F Q({p}, {F}) and that T 
is a convex functional of p(x), as discussed in (a) above. In the course of 
this discussion we describe the behavior of the large deviation functional for 
piecewise constant density profiles. In Section || we calculate optimal profiles 
under various constraints, as discussed in (e) above. In Section we calculate 
the correlations of the fluctuations of the density profile around the most 
probable one and we show that a direct calculation of these correlations 
agrees with what can be calculated from (|2.2| - |2.4| ). Lastly in Section || 
we exhibit a few examples of density profiles for which one can calculate 
explicitly the function F and the value of J-{{p}). 



3 Exact generating function 

For the SSEP with open boundaries as described in the Introduction, the 
probability of a configuration r = {t%, . . . , t^} in the (unique) steady state 



of our model is given by |2d, 21 



Fn[t) ~ WKdTWW) ' { ' 

where the weights wat(t) are given by 

lon(t) = {W\Ylf =l {nD + (1 - n)E)\V) (3.2) 
and the matrices D and E and the vectors \V) and (W\ satisfy 

DE-ED = D + E , (3.3) 

(J3D-SE)\V) = \V), (3.4) 

(W\(aE-~/D) = (W\ . (3.5) 

To obtain the probability of a specified density profile we first calculate 
the sum fijv 1 ,...,iV ra (-Mij M2, ■ ■ ■ , M n ) of the weights ujn of all the configu- 
rations with Mi particles located on the first Ni sites, M2 particles on the 
next N2 sites, etc. The key to obtaining (|2.2| ) is that the following generating 
function can be computed exactly: 

n i Ml 1 • • • ; f^n ) 

l ..Nn 



^ N x \ N n \ 1 n (W\V) 



WW) ' (3 - 6) 
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where the sum is over all Ni,Mi with < Mi < iVj. As shown in Ap- 
pendix |A|, Z can be computed explicitly: 



7 I P0-pl Y +b 
Z = exp 

9 J 



where po and p\ are given by (|T~4|), 

1 



n 

aj^^l - Xi 
L i=l 



(3.7) 



(3. 



7 + a /9 + <5 

and 

g = -pi + poeSLW 1 --^) 

+ ^^L.( e w(i-Ai) _ i )e E 3>i Mi-^). (3 g) 

- ^ A 7 ; 1 



i=l 



Expressions ( |3.7| ) and ( |3.9| ) are the basis of all the calculations leading to 
the large deviation functions. As a first step, note that 

gflj/O-^S^-f W ' P1 V". (3-10) 



<W|V) Vaj-pi-/< 

so that one obtains the normalization factor in (13. II 



_ (W\(D + E) N \V) _ T(a + b + N) 

(w\v) -r(a + b)( P0 - Pl )»- 



4 From the generating function to the large devi- 
ation function 



It is clear from (|3.7|) that Z(X\, ■ ■ ■ X n ; fix, ■ ■ ■ fi n ) is singular on the hyper- 
surface 

g(fix, ■■-,fi n ; Ai, X n ) = 0, (4.1) 



where g is given by (^8|). Let us consider a very large system of N sites 
divided into n boxes of Nx, ■ ■ ■ , N n sites, with PNi,...N n (Mx, ■■■M n ) the proba- 
bility of finding M\ particles in the first box, Mi particles in the second box, 
etc., and let iVj and Mj be proportional to N: iV, = Nyi, Mi = rjiVj = r^iV, 
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i = 1, . . . , n. As explained in Appendix 2, one can show from the definition 
of Z that for large N and fixed yi, ri, 



logP Nl ,... Nn (M 1 ,...M n ) 



N 

n / 



~ iog(p - pi) - Y,yj lo s ~t + r J log a j ' ( 4 - 2 ) 



where the box sizes yj and their particle densities rj are related to the 
parameters px-, ..p n , \±, ...X n by 

dg 



n = rT*kz, («) 

2—ii=\ d log j 
8.9 



, (4-4) 



with all derivatives calculated on the manifold g = 0. 

Equation (|4.2|) gives the large deviation function in a parametric form; 



the 2n + 1 equations (|4.l|) , (|4.3|) , and (|4.4|) determine the 2n parameters 

fix, fj, n , Ai,...,A n in terms of yi,...,y n and ri,...,r n (since yi + h 

y n = 1, the equations (|4.3|) give only n — 1 independent conditions, so 



that the system is not overdetermined) . Note from (3/7) and fl3.9|) that the 
parameters a and b defined by ( |3.8| ) do not appear in the expression for the 
critical manifold and therefore drop out in the large ./V limit, i.e., the large 
deviation functional, like the typical profile, depends only on zq and z\ or 
po and p x . 



4.1 Case of a single box 

One can apply the above results in the case n = 1 of a single box. With 
Ai = A and px = p, equation Q4.1| ) for the critical manifold becomes 

g(p; A) = -px + poe^ 1 ^ + J-^d-A) - 1) = 0, (4.5) 

or more conveniently g(p; A) = 0, where 

Km; A) = Mi -A)- log (^±^1). (4.6) 

V 1 - po + Xpo J 
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The function g may be used in place of g in ([4.3-4.4), because the derivatives 
there are evaluated on the manifold g = 0, so that if the average density 
p = n in the box is given by 



1— pi+Api \ 



+ 



pi A 



Po+Apo / 1-po+Apo 1-pi+Api 



then for large N and M ~ pN, 
1 



iV 



log Pjy(M) ~ — p log A — log 



io & Vl-po+Ap o ; 

(A-l)fa-po) 



Remark: For small fluctuations of the form 

Po + Pi x „ 
p = op <C 1 , 



(4.7) 



(4.8) 



(4.9) 



(|4.8|) gives, again with M ~ piV, 

-ilogP^(M) * " 6 f ? ff 

iV JVV ; 6p + 6pi - 5pg - 2p pi - 5pf 



(4.10) 



For a Bernoulli distribution with average profile (|2.9| ) (the most likely profile, 
as we will see in Section |6|) , one would get 



iio g p„ ( M) * =i«rt! 

3p + 3pi - 2pg _ 2 p oPl - 2pf 



+ 



o ((M 3 ) ; (4.H) 



from this we can compute the fluctuations of M and obtain ( [2.14 ), 



4.2 A finite number n of large boxes 



Let us calculate 



9g 
9 log fj,i 



on the manifold g = 0. From (|3.9| ), we have 



9 log p; 



p oW (l - AOe^^ 1 "^ - p,<^ A 



+p*(l - A,) V — i— (ef^ 1 "^) - l)e^*>i w(1 ~ Ak) . 



(4.12) 
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which becomes, using (4.1) and Q3.9|), 



dg 



d log m 



- Aj) V — L- ( e A«(l-Ai) - l) e E fc>J ^(1-A fc ) + _ (4 13) 

. 1 — Ao 

1>7, J 



One can also show that 
do A,- c 



A; 



SlogAj l-Ai91og^ (Xi-iy 



( e w(l-Ai) _ ijgE^iAyfl-^) .(4.14) 



This parametric form can be simplified by replacing the role of the two 
sequences of parameters Hi an d Aj by a single sequence of parameters G{. 
Let us define the constant C by 



and the sequence Gi by 



c 



G n 



+i 



From the very definition of the GVs it is clear that 

log(Gi/G i+1 ) 



Hi 



1- A* 



and from Q, ( gig ), ( P^D , and ( gig ), we see that 



1- A. 



logCGi/Gf+i) 



^ 1 - Aj C 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



so that both the Hi an d the \ are determined in terms of the sequence Gi, 
the box sizes y, and the constant C. The condition ([4.1|) that 5 = becomes 



Pi . G — Gj + i 



(4.19) 



which gives the extra equation needed to determine the constant C (as 
well as G n+ \). Therefore we are left with n free parameters, the GVs for 
1 < % < n. 
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A more convenient way of writing the Aj's (i.e. ( 4.1§| )) in terms of the 
GiS is 



1 



Vi 



1 



1 



Vj 



Xi-l G i+1 log(Gi/G i+1 ) ru, ^V G i G j+1 J log{Gj/G j+1 



and the condition ( |4.19| ) becomes 



Po ~ Pi 



E Ur. 



{ \Gj G j+1 J logOV^+i) ' 



(4.20) 



(4.21) 



which can be thought as an equation which determines G n+ \ in terms of the 

Once the \ are known through ( 4.2C| ) and fl4,21| ), one gets from (4.4), 
( |4.13| ), and ( |4.14| ) expressions for the rj's and the large deviation function: 



A, 



A; 



Gi — G 



n 



i+1 



l-Xi 
and 

log [P Nl ,... Nn (M 1 ,...M n )} 



1 " A ? ; 



N 



(4.22) 



(4.23) 



i=i 



yi(po - pi) 



Vi log(l - Ai) + yiri log(Aj 



Equations (4.2C-4.23) determine the large deviation function for an any 
specified number n of large boxes, i.e. for fixed yi, ri, and N — > oo. 



4.3 An infinite number of large boxes and the continuous 
limit 

Letting n become large while keeping each box a small fraction of the total 
system, i.e., all the y^s are small or, more formally letting iV — > oo followed 
by n — > oo and yi — > 0, one can introduce a continuous variable x, < x < 1 
and let 



2/1 + 2/2 + ••• + 2/i 



(4.24) 



All the discrete sequences can now be thought of as functions of x with 

Gi = G(xi) ; Xi = X(xi) ; n = p(xi), i = l,...,n. (4.25) 
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Then using extrapolations to make G, A, p smooth functions of x so that 

Gi+i - G t ~ y;G"(x), (4.26) 
one finds that ( f4.20| ) and (4.21) become 



1 



A(x) - 1 G'(x) 



du 



Pi 



r , , , (4.28) 
o G(u) 

At this stage it is convenient to replace the function G(x) by another function 
F(x) defined by 

r x An 

(4.29) 



F(x) = p + f 
Jo 

The expression ( |4.27| ) of X(x) becomes then 

1 F' 2 (x) 



du 



\{x) 



F"(x) 



F(x) . 



(4.30) 



Using the above relations we may rewrite ( 4.22| ) and (4.23) in terms 
of F, obtaining (2.2) and (2.3); the boundary conditions (2.4) come from 
( |4.28| ) . The monotonicity of F follows from the uniqueness of the sign of G 
(see (|4.16|) ; note from (|4.29 ) that G is negative in the case po > p\ that we 
consider). 



5 The large deviation functional in the continuum 
limit 



We derive here some properties of equations (eT3— 2.4 ). We discuss only the 



case 1 > po > p\ > and concentrate on results for piecewise constant den- 
sity profiles, in some instances giving only a sketch of the extension to more 
general p(x). The case in which either pq = 1 or p\ = seems technically 
more difficult; for this case we can show the existence of a solution, but omit 
the proof here. 
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5.1 Uniqueness 



In this section we show that any monotone solution F(x) of ( |2,3D and (|2.4j ) 
is unique. If F and F are distinct solutions then necessarily F'(Q) ^ F'(0), 
since the standard initial value problem for (^), with prescribed values of 
F(0) and -F'(O), has a unique solution. Suppose that F'(0) > F'(0); we 
will show that then F(x) > F(x) for < x < 1, contradicting -F(l) = 
F(l) = p\. For otherwise let y to be the smallest positive number for which 
F(y) = F(y), so that 

F(x) < F(x) , < x < y < 1. (5.1) 

Let = F(x)(l - F(x))/F'(x) and g(x) = F(x)(l - F(x))/F'(x). Then 
(7(0) > g(0) and, from (|2,3|), which can be written in the form g'(x) = 
l-F(x)—p(x), and ([51]), g'(x)-g'(x) = F(x)-F(x) > for < 2 < y, so 
that g(y) > g(y) and hence F'(y) < F'(y), which is inconsistent with (|5.l| ). 



5.2 Piecewise constant profiles and extensions 

In this section we prove existence of a solution F p {x) of ( |2.3| ) and (2.4), 
and show that this solution maximizes G({p},{F}), given by ( |2.5[ ), for a 
piecewise constant density profile 

p(x) = Ti for Xi-i < x < Xi, (5-2) 

where = xq < x\ < ■ ■ ■ x n = 1. We will write y% = Xi — Xi-\. It follows 
from ( p. 3[ ) that F p (x) must satisfy 



i^(s) = >Wi(F,(aO) (5.3) 
for < x < aci, where A{ is a (negative) constant and 

/FY* fl-FV^' 



, r j / V 1 — r, . 

Continuity of F p (x) and -Fp(x) implies that F p (xi), . . . , F p (x n -i) and the 
constants A\, . . . , A n must satisfy 

1 /-^OO t 
Ai = - — - (5.5) 

Vi JF p (xt-i) WW 

and 

^(FpCa^)) = A t+1 i> i+1 (F p (xi)) . (5.6) 
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Conversely, the existence of values F p {x\), . . . , F p (x n -i) and Ai, . . . ,A n sat- 
isfying ( |5.5D and (|5.6|) implies the existence of a solution F p of (|2,3D and 
Q2.4D, obtained by solving (|5.3D on each interval 



So we need to prove that (5.5) and (5.6) can be solved. Now it follows 
from ([5.3D that 



G({p}, {F}) = G ({r}, (F(x ), F( Xl ), F(x n ))), (5.7) 
where {r} = (r%, . . . , r n ) and for {H} = (Hq, . . . , H n ) a sequence satisfying 
Po = H >H 1 >--->H n = p 1 , (5.8) 

we have defined 

1 f Hi dz 



n ( 1 r H i d? \ 

eo(w,m) = E„io g (--/ ffi _ iSM )-io g («- w ). 



(5.8) 



Equation ( p.7| ) suggests that we consider the problem of maximizing Qq. 
Since Qq is continuously differentiable in {H } on the interior of the compact 
domain ( |5.8[ ) and equal to — oo on its boundary, it achieves a maximum at 
some interior point {H*} at which 



dGo 



({r},{H*}) 



1 



1 



for i = 1, . . . , n — 1, where 



A i+1 ^ i+1 (H*) AMH* 



0. 



Ai = - 



dz 



(5.10) 



(5.11) 



Since (|5.10 ) and ( 5.11| ) correspond to ( ^^ ) and (|5.6|) , we obtain a solution 
of the latter equations, and hence a solution F p of ( |2.3[) and ( p.4[ ) by taking 
F p (xi) = H*. Note that the argument above could be used to construct a 



(5.12) 



solution of ( p. 3[ ) and 
da 



from any point at which 
({r},{H}) =0, z = l,...,n-l; 



since the solution is unique (see subsection 5.1), {H*} is the only point 
satisfying ( ^.12|) . 
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It is now easy to verify (2.7) for p(x). For if F(x) is any continuously 
differentiable monotone function with F(0) = po and F(l) = pi, then from 
Jensen's inequality applied on each interval [xj-i, ajj], 



9({p},{F}) < Elftlogf-- / 



r F bi) dz 



log(p - Pi) 



< S u P g ({r},{H}) 

= g ({r}, {H*}) = G({p}, {F p }) = F({p}). (5.13) 

We now discuss briefly the extension of these results to arbitrary profiles. 
A general proof of existence of a solution of ( |2.3j ) and (|2.4| ) may be given 
which is independent of the arguments above: if p(x) is continuous then 



Theorem XII. 5.1 of [35] implies immediately the existence of a solution F p of 
( |2.3|) a nd (2.4), and for measurable p(x) only slight modifications of the proof 
in [^g] are necessary. The uniqueness theorem for the initial value problem 
for ( |2.3| ) implies that the derivative of the solution must be everywhere 
nonzero, so that the solution is monotonic. Equation ( |2.7D may be verified 
for arbitrary p{x) by a limiting argument from the same result (proved 
above) for piecewise constant densities; the key idea is to show that F p and 
Fp are, in the uniform norm, continuous functions of p in the L 1 norm. 

5.3 Convexity of ^({p}) 

Finally we show that J-{{p}) is convex. Recall that F p is the solution of 



(2.3) and (2.4) corresponding to the profile p. From (2.6) 

5G 



5F{x) 



{p},{F P } 



0, 



(5.14) 



we have 

5 2 g 



5p(y)5F(x) 
and therefore 
8 2 F 
Sp(x)Sp(y) 



+ 



5 2 Q 



{p},{F P } Jo 5F(u)5F(x) 



6W ({p})du = 0, (5.15) 



(W) 
5 2 G 



S ( 8Q 
5p(x) \8p(y) 



{pUFp}, 



5p(x)5p{y) 

fl rl 



{p},{Fp} 

5 2 G 



(5.16) 



o Jo 5F(u)5F(w) 



{p},{F P } $p(x 



6Fp{U \{p}) 6 44^({p})dnd W . 



Sp{y) 
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Since G({p}, {F}) has a maximum at F p the second term in this expression 
is positive semidefinite. The first term is positive definite: 



5 2 Q 
5p(x)5p{y) 



5(x - y) 

{p},{F p } p(x)(l - p(x))' 



(5.17) 



6 Optimal profiles 

The convexity of J~({p}) established above implies the existence of a unique 
global minimum, corresponding to the optimal or most likely profile p(x). 
From (p.8|) it follows that p(x) must satisfy 



log 



p(x)(l-F(x)) 



[(1 - p(x))F(x) 



0, 



(6.1) 



where F(x) is the solution of ( jO| ) and (|2.4D corresponding to p{x). Equation 
( |6.1D leads immediately to 

F(x) = p(x), (6.2) 

and with 
then yield 



this implies that F"(x) = 0; the boundary conditions (|2.4| ) 



F(x) = p{x) = p + (pi - p )x, (6.3) 

verifying that the optimal profile is as given in (|2.9| ). 

We may also ask for the most likely profile p(x) under a constraint of 
the form 



ip(x)p(x) dx = K, 



(6.4) 



with fixed weighting function ^(x) and constant K. From ( |2.8| ) one sees 
that p(x) must satisfy 



p(x){l-F(x)) 
(l-p(x))F(x) 



exp[cip(x)} 



(6.5) 



for some constant c, where again F(x) is the solution of (2.3) and (J2T 
corresponding to p(x). 

In particular, we may impose the constraint of a fixed mean density p* 
by taking ip(x) = 1, K = p*. Then from (|6 . 5|) , 



p(x)(l-F(x)) 
(l-p(x))F(x) 



(6.6) 
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and we find from fl2.3| ) that F must satisfy 

F"(x) e c - 1 



F' 2 (x) 1 + (e c - l)F(x) ' 



(6.7) 



The solutions of (6/7) have the form F(x) = A + Sexp(— 6x), which with 
the boundary conditions fl2.4|) leads to 



F(x) = (Pi-W^-fa-f*)*- * . m 

from fl2.3|) one then finds an exponential profile 

(jc) = [(1 - P0 )e- e - (1 - Pl )][( Pl - poe- 9 )e 9 * - ( Pl - p )] 

(1 - e~ e )(pi ~ Po) 

Here is a free constant which must be chosen to satisfy the constraint 
of mean density p*. If p* = (po + pi)/2, the expected (and typical) total 
density in the stationary state, then = and ( |6.9| ) reproduces the optimal 
profile p(x). 

Remark: (a) Unless = or po = 1, 

p(0) = limp(x)jtp . (6.10) 

cc\0 

Thus when we constrain the total number of particles to be different from 
its most likely value, the optimal profile is discontinuous at x = (unless 
po = 1). Similar conclusions hold at x = 1, where a discontinuity occurs 
unless = or p\ = 0, and, by symmetry, in the case po < p\. 

(b) In the limit p\ /* po, one should scale ~ po — p\ in order to obtain 
a meaningful answer. The resulting optimal profile is constant, with an 
arbitrary density po + po(l - po) / (po - Pi) ■ 

(c) If is chosen such that 

1 - Po + PoA 

one can recover 

and Q; (fOj) by calculating $p(x)dx from Q, 
and dJ) from Q, Q and ( |OD . 

Finally, we may also impose simultaneously several constraints of the 
form (|6.4j) . For example, we can decompose the system as the union of 
disjoint boxes and then fix the mean density in some of these, imposing 
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no constraints in the remainder. Then ( |6.1| ) will hold in the unconstrained 
boxes, so that F(x) = p(x) will be linear there, and ( |6.6| ) will hold in the 
constrained boxes (with a box-dependent constant c), so that F and p will be 
exponential there; the specified mean densities in the constrained boxes and 
the requirement of continuity of F and F' at the box boundaries then com- 
pletely determines F and hence p, which will, in general, be discontinuous 
at the boundaries. 

Suppose, for example, that we require that the density vanish for < 
x < xo, with no additional constraint. Then 

n*) = { 1_( n" Po) C^o. n l X ~ X(h (6-12) 
\ pi - (1 - po)Be Bx °(x - 1), tfx>x , 

and continuity of F at x = xq requires that B satisfy 

l^± = {l + B(l-x )}e Bx °. (6.13) 
1 — Po 

From Q2,2| ) we find that 3~({p}), which in this case is simply the probability 
that p(x) vanish for < x < xq, is given by 

H{ P }) = - log (^) + (1 - *.) log ( 1 + ^% ) ■ (6-14) 



If now we take po — ► 1 with p\ and xq fixed, ( 6.1 3j ) implies that B ~ 
— Xq 1 log(l — po) and therefore 

T({p}) ~ x log(- log(l - po)) • (6.15) 

Thus the probability of finding zero density in the box x < xq when po = 1 
is super-exponentially small. 

7 Small fluctuations in the profile 



We have already given in (4.10) the formula for the probability of small 
fluctuations in the total number of particles in the system. In this section 
we compute, directly from the large deviation functional 3~({p}), the co- 
variance of small fluctuations in the profile around the stationary profile 
p(x) = (1 — x)po + xp\ (|2.9|). We show also that the result agrees with 
with a direct computation of the correlation functions from the microscopic 
measure describing the SNS (note that the possibility that these two com- 
putations might disagree had been conjectured in 0]). 
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7.1 Small fluctuations from large deviations 

For a small fluctuation of the density profile around its optimum p(x) , 

p(x) = p(x) + e(x) , e(x) « 1 , (7.1) 
there will be a corresponding variation of F(x), 

F(x) = p(x) + <f>{x). (7.2) 



From ( |2.4|) we have 0(0) = (f)(1) = 0, and from ( |2.3| ) </>(x) will satisfy, to first 
order in e(x), 

e(x) = 0(a? + -r^2 <P (x). (7.3) 



From (|2.2j ), again to lowest (quadratic) order in e(x), 

HW) = -lJ^{^- ^^ t"(*f}. (7.4) 
One can rewrite ( |7.3| ) as 

-■I 



o 



e(x)= / dy C(x,y) </>"(y) , (7.5) 

with C(x,y) given by 



C(x,y) = -(l- x )y8(x-y)- X (l-y)6(y-x) + p{x){1 _ /2 P{x)) 5(y-x), (7.6) 
where 6(x) is the Heaviside function, or with a more compact notation, 

C(x,y) = A-\x,y) + p{x){l ~ p{x)) 5{x - y). (7.7) 

piz 

Here A is the second derivative operator with Dirichlet boundary conditions 
at x = and x = 1. The expression Q7.4|) can then be written as 



^({P» = T^o I dx I dycf>"(x)C(x,y)ct>"(y) 



2p' 2 Jo 
1 



2p' 2 Jo 



/ e(x) C7" A (x,i/) e(y) , (7.8) 



where we have used (7.5). This formula expresses the fact that, for large 
N, density fluctuations are approximately Gaussian with covariance matrix 
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p' 2 C(x,y)/N. Using the explicit expression Q7.6Q of C, this leads to the 
formula for the correlations (e(x)e(y)): 

(e(x)e(y)) = 



N 

^{-p ,2 [(l - x)y9(x - y) + (1 - y)x% - x)] 

+ p»(l-p(xMx-!,)}. (7.9) 



This result was obtained earlier by Spohn || , and shown there to agree with 
the results from fluctuating hydrodynamics. 

7.2 Microscopic correlation functions 



One may also compute correlations directly from from the algebra (0-3.5) 



and the normalization factor ( p. 11 ). By recursions over i and j one finds for 
1 < % < N that 

(W\(D + E) i ~ 1 D(D + Ef-^V) 

= Po (W\(D + E) N \V) -(a + i- 1){W\(D + E) N ~ l \V) , (7.10) 

and for 1 < j < i < N that 

(W\(D + Ey~ x D{D + E) i ^- 1 D{D + E) N -*\V) 
= pl(W\(D + E) N \V) 

-p (2a + i + j- 2)(W\(D + E) N ~ l \V) 

+(a + j-l)(a + i- 2)(W\(D + E) N ~ 2 \V) . (7.11) 



Using ( |3.11 ), one obtains from ( 7.10| ) the average profile 



a + b + JS — 1 

and from ( [7.1 1| ) the truncated (microscopic) correlation function: for j < i, 

i \ i \i \ i ^2 (a + j-l)(b + N-i) 

{rjn) - Win) = -0* " po) {a + b + N _ ma + b + N _ 2) ■ ( 7 - 13 ) 



Equations §Ta% and (ffl3|) agree with the corresponding formulas in 



where for rates satisfying a + (3 = a = j + 5 = b= 1, random walk methods 
were used to calculate the one particle and pair correlation for this system. 
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In the large N limit, ( [7.12 ) yields the linear profile (2.9), and by summing 
( |7.13| ) over i and j, one recovers ( [2. 14 ) Moreover, we may recover (|7.9|) : if 
we divide the system into boxes of size Ax (with 1 <C Ax^ 1 <C N) and 
write, for x and y multiples of Ax, 



e(x) 



N(x+Ax) N(y+Ax) 

n-p(x), <y^-j^ ^ T i-p(y)' ( 7 - 14 ) 



NAx 



i=Nx 



i=Ny 



then 



(e(x)e(y))-(e(x))(e(y)) 



1 



iV 2 Ax 2 



N(x+Ax) 

5 x , y £ ((n)-(Tj) 2 ) 

i=Nx 



N(x+Ax) N(y+Ax) 

+ E E (fonJ-foXTi))}, (7.15) 

j=Nx i=Ny 



and this, with (|7.12|) and (OS), yields O). 



8 Examples 

Although in general the calculation of the large deviation function J-({p}) 
for a given profile p(x) is not easy, since it requires the solution of the non- 
linear second order differential equation ( |2.3| ), there are nevertheless some 
cases for which it can be done. Piecewise constant profiles were discussed in 
Section ||.2; in particular, for a profile of constant density r the large devia- 
tion functional is given by ( |2.12| ), and this function can be found explicitly 
when r = or 1 (in which case F(x) is exponential) or r = 1/2 (when F(x) 
is sinusoidal). F({p\) can also be found explicitly in some cases for the 
optimal profiles under constraint, discussed in Section |6| 

One may also construct examples by specifying F(x) and obtaining 
p(x) from (|2.3D and !F({p}) from One must check in each case that 

< p(x) < 1, since this is not guaranteed by (p.3|). In the remainder 
of this section we give two examples of this type which address the ques- 
tion: for what profiles p(x) is ^({p}) infinite? If 1 > po > p\ > then 



3~({p}) < 00 i this follows from (2.2) and the fact that, by the uniqueness 
theorem for the initial value problem for ( |2.3[) (see remark at the end of 
Section 5.2), F'(x) cannot vanish. Suppose then that po = 1 (analysis of the 
case p\ = is similar). The examples below suggest that ^({p}) = oo when 
lim^^o p{ x ) = and this limit is approached faster than any power of x. 
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This is also supported by the example at the end of Section ||, which shows 
that ^({p}) = oo if po = 1 and p(x) vanishes identically on an arbitrarily 
small interval < x < xq. However, we have not formulated any sharp 
conjecture. 

In both examples we take po = 1- F° r the first, define 

F(x) = l-{l- Pl )e c( - l - l l xn \ (8.1) 

with c a positive constant. Note that F is monotone decreasing and satisfies 
the boundary conditions fl2.4|). Then 



n c 

F\x) = - — (l-F(x)), (8.2) 



F"i:r) - ( ^-V f 1 _ ^±l x » ] (1 - FCx)), (8.3) 



and hence from (127 



x n+l ii cn 



n 4- 1 

p( x ) = !iXi^F(x). (8.4) 

cn 

Thus p(x) vanishes like x n at x = 0, while lim ;c _ + i p(x) = pi(n + l)/cn. 
Certainly we may choose c (e.g., c = (n + l)/n) to ensure that < p(x) < 1 
for all x (note that by choosing c sufficiently small we see that ( |2.3|) for an 
admissible F does not guarantee < p(x) < 1). If we then write ( |2.2[ ) in 
the form 



t(i /Mm, p(x)(i-F(x)) 

F({p» = y o dx|p(x)io g (1 _ p(x))F(x) +iog 



p(x)(l-F(x)) , , (l-p(x)) 
(1 - Pi) 



+ log (T^ ^ 



then (§3]), (|OD , and fyj) show that .F({p}) < oo. 
Finally, we again take po = 1 and c > and define 



Then 



F(x) = !-(!- pi)e c(e - elA) . (8.6) 



F\x) = -¥—(l-F{x)), (8.7) 

x z 



:r 2 



26 



and 

p( x ) = t±—± e -y x F(x). (8.9) 
c 

Now p{x) vanishes faster than any power of x at x = 0. Again we may 
choose c (e.g., c = 3) so that < p{x) < 1 for all x; if we then write (2.2) as 



in ( p.q ) we see from (|8.6| ), (|8.7| ), and (|8.9| ) that although the first two terms 
are finite, the last contributes Jq 1 dx/x, so that ^({p}) = oo. 

9 Conclusion 

In the present work, we have seen that the large deviation functional ^{{p}) 
of the density profile p(x) can be calculated for the SSEP, starting from the 
known weights of the microscopic configurations. We find that the large 
deviation function is in general finite, although there are counterexamples 
(see fl2.12|) and section ^) . 

The simplest way we found to write our results is a parametric form 
(f2~2 -2.4) in which both the large deviation function and the density profile 



are expressed in terms of a monotonic function F(x) varying between pq and 
pi, the densities of the two reservoirs at the two ends of the system. 

An interesting question posed by the present work is how the additiv- 
ity property of the free energy of equilibrium systems is modified in the 
nonequilibrium case. We first note that if a system of N sites, in con- 
tact with left and right reservoirs at chemical potentials corresponding to 
densities p a and p , is described by a macroscopic coordinate x satisfying 
a < x < b (rather than < x < 1), then the probability P of observing a 
profile p(x), a < x < b, satisfies 

N 

logP~-- F[a,b}({p};Pa,pb), (9.1) 



where 



F[a,b]({p}'-> Pa, Pb) = J dx!^p(x)log (yr^y) 



+(1 _ pw)Iog( l^) +log ((^« u , ( , 2) 



and F{x) is related to p{x) on the interval a < x < b by (|2.3| ) with the 
boundary conditions 

F(a) = p a , F(b) = Pb . (9.3) 
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(Note that F'(x) has the same sign as Pb — Pa, so that the argument of the 
log is positive in the last term of (|9.2[) .) 

Now consider two systems, of Nu and N(l — u) sites respectively (with 
< u < 1), with the left system in contact with left and right reser- 
voirs at chemical potentials corresponding to densities po and p m and the 
right system with reservoirs corresponding to densities p m and p\ (with no 
other relation between them, in particular with no particles jumping di- 
rectly from one system to the other). If p{x) is a profile on < x < 1 
and p^(x) and p^ 2 \x) are its restrictions to the intervals < x < u 
and u < x < 1, then we would like to relate the probability P/v of ob- 
serving p(x) to the probabilities P$ and pjp of observing p^\x) and 
p( 2 \x) in the two subsystems. If T were truly local we would have P/v — 
pffpffi. The most naive generalization of this idea here would be that P/v — 
sw P Pm ( P N P N ) or, since from logP^ ~ -NF [0tU] ({p (1) }; Po, Pm) and 
logP^ 2) ~ -iV^[ U) i]({p( 2 )};p OT ,pi), equivalents that .P[o,i]({p};po, Pi) = 
mf p™.[ jr [o,n]({p (1) };Po,Pm) + F[u,i]{{p {2) Y, Pm, Pi)]- One can check, however, 
that this is not true. 

Instead, we observe that if we define a "modified free energy" 7i by 

{{P}; Pa,Pb) = F[a,b] ({p};Pa,Pb) + Q> - o) log ' ( 9 " 4 ) 

then we obtain from ( p. 7] ) the "additivity" property 

W[ ,1]({P}; PO, Pi) = SUp{W [0iU ]({p (1) }; /3 , Pm) + H[ U ,1] ({P (2) }; Pm, Pi)}- 

(9.5) 

The occurrence in ([O]) of the supremum rather than infimum is a conse- 
quence of ( |2.2| ), but we do not understand its physical basis at this time. 

It is perhaps surprising that the additivity property ( |9.5| ) suffices to 
determine the large deviation functional completely. For suppose we divide 
our system as above, but into n rather than 2 subsystems, with division 
points = xo < x\ < ■ ■ ■ < x n = 1, and denote the corresponding reservoir 
densities by po = Po > Pi > • ■ ■ > P n -i > F n = p\. Then iterating (|9.5| ) we 
find that 

•^[o,i]({p};Po,Pi) 

= sup E I T {x . _ 1>Xj] ({p}; Fj-x, Fj) + yj log ( p' 1 ~ f J ) 1 ,(9.6) 

fi.-A-!^ 1 ' J \(po-pi)yjJ J 
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where yj = Xj — Xj—i. When n is large and all yj are small, 



F i -F j - l ~y i F'(x j ) 



(9.7) 



moreover, since Fj 



that is, the two reservoirs at the ends of the in- 



terval \ ] have essentially the same chemical potentials, we can expect 

each this subsystem to be in equilibrium, so that from ( 2.10f) , 



p(Xj)\o 




+ (l-p(a?i))log 




Vi 



By substituting fl9.7| ) and Q9.8Q into ( |9.6| ) and taking the limit n 
yj — ► for all j, we obtain our basic result (fE?|). 

Clearly it would be of great interest to give a direct derivation of the 
additivity property ( |9.5[) , and to know if this property is limited to the 
SSEP or whether similar relations hold for more general systems. 
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A Derivation of ( |57T1) and ( CT ) 

A.l A first consequence of ( |3.3|) 

Let us first prove that if D and E are operators satisfying ( p.3[ ), 

DE-ED = D + E, (A.l) 

then 

c xd+ve = ( ( x - y) eV \ E f ( x ~ y) eX \ D _ (A2 ) 

\xeV — ye x J \xeV — ye x J 
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Equation ( |A.2| ) and similar equations below are to be interpreted in terms 
of formal power series in x and y. 



One can easily check from (A.l) that for all p > 

DE P = {E + IfD + E(E + If - E p+1 , (A.3) 
which means that for "arbitrary functions" f(E), one has 

Df(E) = f(E + l)D + E[f(E + 1) - /(£)] . (A.4) 
Then if one tries to write e z ( xD +v E ) under the form 

e z(xD + yE) = e tE e uD ^ ^ 

one gets that 

^ = x(e t -l) + y, (A.6) 

- = xe\ (A.7) 
and by integrating over z, one obtains ( |A.2| ). 

A. 2 Other consequences of ( |3.3| ) 

Using ( |A.2| ), one derives easily the following three identities: 

eXD+ yE = f xe x -yey \ D f xe*-yey \ E 
\{x-y)ey) \{x-y)e x 



Then combining ( A.8 ) and ( |A.S| ), one can also show that 

(po - Pi)e x ^ poE-(i-po)d 



e xE e yD 



(A. 



= ( t Y ( t Y (A 9) 

\ e x + e y _ e x+y J \ e x +e y_ e x+yJ ' v. ■) 

W-(*±£±)° (*±*zl)' . (A ,„) 



1 - (1 - p )e x - Px ev 

(ft)-px)c y n 
1 - (1 - p )e x - piey 
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, (A.ll) 



and that 



e uE g vD e xD+yE 



(x - y)e u+y 



(x - y)ey + y(e v +y - e v+x ) J 
j (x — y)e v+x 



D 



(x - y)ey + y(e v+ y - e v+x ) 



.(A.12) 



A. 3 Derivation of ( P77|) 

As a consequence of ( A.12f) we see that if {x n } and {y n } are two sequences 
and if the sequences {u n } and {v n } are defined by the recursion 



e M„+i-E giin+iD _ e u n E e v n D e x n+ iD+y n+1 E 



with the initial condition 

u = v = 0, 

then the general expressions of u n and v n are given by 



?£,7=iM- x * + — - — (e Vl ~ Xl - \ )(^ ■ ' 



-i -l 



(A.13) 
(A.14) 

(A.15) 



En 



n 

E" , w + ^ _ ^ ( e »-^ _ i) e E 3>1 v,-^ 



i=l 



Therefore 
where u n and v n are given by 



e AtiAiD+^iE _ _ _ e ^ k X k D+ii k E _ e u n E ^nD 



(A.16) 
(A.17) 



n -i 
• -, A,; — 1 



i=l 



(A.18) 



„/<i(l-Aj) 



l)e 



E,- A • 



(A.19) 
n -1 



e Er=i 

n 1 

i=l Ai " 1 

Then using ( A.ll| ), one can rewrite ( A.17| ) as 

gMiAiD+MiE , , , el i k \ k D+v k E = e U n [-D+p (D+E)] e V n [D- Pl (D+E)\ ^ (A.20) 
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where 

e V n = Pl-PO ^ e Un = e K l +X;r=l^( 1 - A »), (A.21) 

9 

with g given by (pT9). Lastly, with po, Pi, a, and 6 given by (1.4) and (j3] 



the algebraic rules ( 3~4|) and (|3.5| ) can be written as 



[D-pi(D + ^)]|V*) =6|V) , (A.22) 
(W\[-D + p (D + E)} =a{W\ , (A.23) 



and one obtains Q3.7| ) from ( A.2C| ). 



B Derivation of (jOH^) 

Consider a system of A r sites, divided into n boxes, with a fugacity Ai on 
the N\ first sites on the left, then A2 on the next N2 sites and so on, A n on 
the last N n sites. Clearly we have 

N = N 1 + N 2 + ... + N n . (B.l) 

Let us define fi by 

(W\(XiD + E) Nl ...(X n D + E) N " \V) 



(W\V) 



\Mi x M n ^ l N 1 ,...,N n (M 1 ,. . .,M nJ 
" 2. A l ■•• A n » ( B - 2 ) 

(see (|3.6|) ), and f2o as in ( |3.11|) , 

= (R3) 

Here we have suppressed the dependence of Qq on N and of f2 on Ni, . . . , N n 
and Ai, . . . , X n . Suppose that for large Ni, .., N n , the quantity defined by 
( |B.2| ) has the following behavior 

Q ~ e NHVim,---yn;Xl,..Xn) Nl l ^ Nn ] j ( R4 ) 

where 

K = f • (B.5) 
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Clearly one has that the average density T{ of particles in box i, in the 
ensemble with fugacities Ai, . . . , A n , is given by 

tVl ~ N " iVSlogA, - SlogAi • 1 bj 

If we assume that the distribution of the Mj is strongly peaked near their 
mean values then 

n = n ]T Af 1 ---A* / "P JVl ,..., iVn (M 1 ,...M n ) 

M u ...,M n 

« n X[^ N ■ ■ ■ \ r n ^ N P Nl _ Nn (r iyi N, ...r n y n N) , (B.7) 



where P/Vi,..jv n (-Mi, ...M n ) denotes the probability computed in the ensemble 

(B.8) 



in which Aj = 1 for all i. From (B.4) and 
logP Nlt ... Nn a(Mi,...M n ) 



N 

~ h(y 1} ...y„; Ai, ...A„) - nyi log Ai... - r n y n log \ n + K , 
with Mi = riyi and K given by (see (|3,11|) ) 

K = + 1 ^ - -pi) + ± m log y, . (B.9) 

i=l i=l 

We thus obtain Pjv 1) ..jv„(Mi, ...M n ) in a parametric form (|B.g| ): as we vary 
the Aj's, the densities r{ in the boxes vary according to ( |B.6| ). 

Let us now see how we can extract the function h which appears in ( |B.4| ) 
from (|3.6|-|3.9D . If we consider the generating function (|3.6|) , 



V y Pi 1 ^ (W\(X 1 D + E) N K..(X n D + E)^\ V ) 

^ "' ^ Nir'NJ. (W\V) ' { } 

N 1= N n =0 1 n V v \ v I 



and use the asymptotic form ( |B.4| ), we see that Z becomes singular along a 
manifold 

9 = 0, (B.ll) 

where 

g= max [h(yi, ...y„; Ai, ...A n )+yi log/ii + ...y n log//„]. (B.12) 

j/l,...y„with y 1 +...y n =l 
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Along this manifold, we have obviously 



h 



(B.13) 



so that ( |B.8| ) and ( |B.9D reduce to ( |4.2D . Moreover we see from ( gig ) that 

1 * (B.14) 



Cdlog^ ' 



where the constant C is a Lagrange multiplier associated to the constraint 
that J2i Vi = 1) an d this establishes (4.3). Las tly fl4.4p follows simply from 
the fact that along the hypersurface given by ([B.11D and (]Bl2|) , one has 



Oh 



dlogXi dlogXi ' 
which shows that ( |4.4| ) follows from ( |B.6| ). 



(B.15) 
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